Extended Impurity Potential in a d x 2_ y 2 Superconductor 
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We investigate the role of a finite potential range of a nonmagnetic impurity for the local density 
of states in a d x 2_ y 2 superconductor. Impurity induced subgap resonances are modified by the 
appearance of further scattering channels beyond the s-wave scattering limit. The structure of the 
local density of states (DOS) in the vicinity of the impurity is significantly enhanced and therefore 
improves the possibility for observing the characteristic anisotropic spatial modulation of the local 
DOS in a d x 2_ y a superconductor by scanning tunneling microscopy. 
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A lot of attention has been focussed recently on the role 
of impurities in unconventional superconductors. The 
reason for this interest is that impurities modify the su- 
perconducting properties in a way which is characteristic 
for the pairing state and thereby serve as a diagnostic 
tool for its identification. In particular, magnetic impu- 
rities act as strong pair breakers for superconductors in 
a spin singlet pairing state while already non-magnetic 
impurities are pair breakers in superconductors with a 
nontrivial phase of the pairing amplitude. The latter 
kind of defects or impurities have only little effect on 
the transition temperature and the superfluid density in 
conventional s-wave superconductors as understood from 
Anderson's theorem jjj. 

Regarding high-T c superconducting materials evidence 
has accumulated for an anisotropic energy gap most likely 
of d x 2_ y 2 symmetry 0. For this pairing state nonmag- 
netic impurities produce a finite lifetime for quasiparti- 
cles near the nodes of the gap and a finite DOS at low 
energies [§|-[§. Indeed, the measured low temperature 
properties of high-T c materials show a remarkable sen- 
sitivity for impurity effects. Examples are the T 2 vari- 
ation of the low temperature penetration depth |q] and 
the Knight shift in Zn doped YBCO 0, and the u? to w 
crossover in the low frequency B\ g Raman intensity 

A general feature of pair breaking impurities in s-wave 
superconductors is that they lead to bound quasiparticle 
states in the energy gap Q . Similarly, for a <i-wave super- 
conductor with a particle hole continuum and a DOS ex- 
tending linearly to zero frequency, low energy resonances 
are created with a highly anisotropic structure [jio] [L2|. 
Increasing the impurity concentration these subgap reso- 
nances form bands which give rise to the afore mentioned 
finite DOS at zero energy. The local structure of the 
resonant states in the vicinity of the impurity is char- 
acteristic of the pairing state of the superconductor. It 
has therefore been suggested that probing these impurity 
states by scanning tunneling microscopy (STM) offers a 
direct way for measuring the anisotropic structure of the 
pairing state |13|. However, the predicted effect may be 



too small to detect with the current STM spatial resolu- 
tion and thus any advances that could enhance a charac- 
teristic signal would be helpful towards making STM a 
feasable probe to detect the gap symmetry fll4j . 

Systematic studies in high-T c materials have been per- 
formed by substituting Zn on the planar Cu sites and 
thereby suppressing the local moment. Model calcula- 
tions, as initiated by the early work of Annett et al. and 
Hirschfeld and Goldenfeld ||, which assume a d x 2_ y 2 gap 
function with unitary s-wave impurity scattering have 
provided a consistent explanation for the experimentally 
observed transport properties. The origin of the unitary 
nature of this scattering from local defects has been as- 
cribed to strong electronic correlations. While this issue 
is unsettled, it has in fact been suggested that a Zn im- 
purity may not simply act as a local moment vacancy 
because it induces a magnetic moment in its surrounding 
in the C11O2 planes Jl5| | . 

Transport data on Zn doped metallic YBCO have fur- 
thermore been used to infer a spatially extended nature 
of the effective impurity potential. In particular, a scat- 
tering cross section diameter larger than the Cu-Cu dis- 
tance has been deduced from the residual resistivity jl6) . 
Theoretically, for a model calculation of static vacancies 
in a Heisenberg antiferromagnet it has been shown early 
on that the ordered staggered magnetic moment is en- 
hanced within a few lattice spacings around the vacancy 
p7[ . Therefore, even in the hole doped system when spin 
correlations are dynamic and have a finite range the spa- 
tially extended structure of impurity effects is expected 
to persist. Electronic correlation effects thus provide a 
natural origin for the transport data. 

A series of model studies has been performed in re- 
cent years on the properties of dirty cZ-wave supercon- 
ductors. Most of the analytical studies have assumed 
s-wave impurity scattering while numerical simulations 
of impurities in Hubbard or t—J type models have in fact 
demonstrated the expected effective, dynamically gener- 
ated finite impurity potential range ji8[ . As a conse- 
quence higher order partial wave scattering channels be- 
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FIG. 1. Model for a nonmagnetic extended impurity on a 
square lattice. Vb is the potential at the impurity center and 
Vi the potential at the neighboring (nn) sites; ti is the change 
of the nn hopping amplitude of the host value —t. 

come important and lead to resonances in each channel. 
Impurity models for cuprate superconductors thus need 
to take into account the extended nature of the scattering 
centers |Q. 

In this letter we explore the T-matrix in a mini- 
mal, physically transparent extended impurity potential 
model in a d-wave superconductor. In a partial wave 
decomposition the single impurity problem is solved ex- 
actly. Subgap resonances in the local DOS are strongly 
enhanced by a finite range of the potential and may thus 
be the origin for the unitarity limit scattering which was 
favoured in previous dirty <i-wave studies for cuprate su- 
perconductors. This enhancement of the resonances im- 
proves the possibility for the detection of their highly 
anisotropic spatial structure in STM spectroscopy. Also 
the impurity induced zero energy DOS is found to be 
significantly increased for finite impurity concentrations 
as compared to the s-wave scattering limit for a local 
impurity potential. 

Specifically we model the scattering from impurity cen- 
ters located at the sites {1} on a square lattice 
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(see Fig.l). Here, S connects to the nearest neighbor- 
sites, Vo and V\ < Vq are potential strengths at the im- 
purity center and its nearest neighbor sites (nn), respec- 
tively, and ti is an additional contribution to the nn hop- 
ping amplitude which we assume to be of opposite sign to 
the pure host value —t. By focussing on a two parameter 
model for the specific parameter relation V\ — a 2 Vo/A 
and ti — aVo /A in Eq. ([[]) Fourier transformation leads to 
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where V" k ^ = 1 + cve lk '^. The virtue of the parame- 
ter choice becomes apparent in the factorization of the 



scattering potential which allows for an algebraic solu- 
tion of the single impurity T-matrix. a is the control 
parameter for the extension of the potential to the nn 
sites. The single particle Green function in the pres- 
ence of the impurity is G(k, k',w) = Jk,k'G°(k, k', u>) + 
G (k,w)T k!k '(cj)G°(k',w). All quantities are 2x2 ma- 
trices in Nambu particle-hole spinor space. The clean 
Green function is given by [G°(k, = uitq — A k fi — 

6c t 3 = Ei=o,i,3^ff«] _1 ' wn ere f; (i = 1,2,3) are the 
Pauli matrices, frj is the unit matrix, and £ k is the quasi- 
particle tight binding energy relative to the chemical po- 
tential. Specifically, we consider a superconductor with 
a d x 2_ y 2 gap function A k = A(cosk x — cosk y )/2. 

In order to explore local properties in the vicinity of 
an impurity we consider first the problem of a single im- 
purity located at the origin 1 = 0. The momentum de- 
pendent T-matrix is obtained from the set of equations: 
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where underlined quantities are 4x4 matrices with re- 
spect to the nn site coordinates next to the impurity 
center; note that their matrix entries contain noncom- 
muting Pauli matrices and therefore the matrix inversion 
in Eq.(4) is nontrivial. 

Using the square lattice basis functions = 
i (cosk x ± cosfcj,) and 'j(p 1 'p 2 ') = I (sinfe^ ± sinfc^) and 
thus decomposing into s— , p—, and d-wave scattering 
channels we obtain the algebraic result for the T-matrix 
in the form 

Tk, k /M = y U k D(a;)U k ^3 , 

U k = [(1 + art) , iajl 1 , ia 7 £ 2 , a^] (6) 

where the matrix D,(u>) is given by 
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The structure of the matrix D_ apparently implies a mix- 
ing of the s- and d- and the two p-wave scattering chan- 
nels, respectively. The matrix elements in Eq.(7) are de- 
termined by the 2x2 matrix functions 
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FIG. 2. Frequency dependence of the local DOS at r = 
(1, 1) at Vo — St for an extended (a = 0.4, bold solid line) and 
a local (a = 0, dashed line) impurity in a d x 2_ y 2 supercon- 
ductor. The impurity center is at the origin r = (0,0). The 
thin solid line shows the pure DOS. 
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FIG. 3. Distance dependence of the local DOS at the 
resonance frequency uj = +ui re3 along the diagonal (main fig- 
ure) and horizontal (insert) direction away from the impurity 
center at r = (0, 0). 
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momentum dependent coefficients 

Sk = 4 + 2cv 2 + 807^ + 4a 2 cos k x cos k y — p^. 
c?k = 2a 2 — 8a 2 cos k x cos k y — p^ , 
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with the 
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In the limit a — > the T-matrix becomes momentum 
independent and reduces to the known s-wave scattering 
result of a local impurity potential po|| . Note, however, 
that for finite a the T-matrix T — y\._ Q Tjfi has compo- 
nents in fo and all three Pauli matrices while T\ = = 
for a = 0. Given the explicit algebraic solution of 
Tk,k'(w) the local DOS then follows from 
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r G 11 (k,k', w )e ik ' r (13) 



with u> + i0 + implied. 

We have evaluated the local DOS using a free tight 
binding band £k = — 4< 7 £ + At' cos k x cos k y — fi, with 
t = 1 as the energy unit and t' = 0.4, for a band filling 
n = 0.86 and a d-wave gap amplitude A = 0.3. Yet, the 



conclusions are insensitive to changes in this parameter 
choice. The potential at the impurity center is chosen as 
twice the bandwidth Vo = 8. All results for the extended 
potential, i.e. finite a, are compared to the local potential 
case a = 0. 

In Fig. 2 we have plotted the frequency dependence of 
the local DOS at a next nearest neighbor (nnn) site of 
the impurity center showing a sharp resonance at a low 
subgap frequency uo res . The p-wave scattering channels 
are too weak to contribute an additional structure and 
the single resonance has a mixed origin in the s- and 
d-wave channels. Quite remarkably, at this nnn site the 
strength of the resonance is enhanced by a factor 4 when 
a is increased from to 0.4. Comparably weaker reso- 
nances at ±iiV es with weaker a induced changes appear 
along the (1,0) and (0,1) direction and this anisotropy 
directly reflects the nodal structure of the d x 2_ y 2 gap 
function [flC 



In Fig. 3 we s 



13 



ihow the oscillatory distance dependence of 
the local DOS along the diagonal (and, in the inset, hor- 
izontal) direction at the resonance frequency 
The enhancement due to a finite range of the potential 
is largest near the impurity center where the resonance 
is strongest. In comparing the results for the diagonal 
and the horizontal direction we recognize that the direc- 
tion dependence of the resonant structure is amplified by 
the finite potential range thereby improving the possibil- 
ity for detecting the anisotropic local DOS structure by 
STM in the vicinity of the impurity. 

While this is already a very promising result for lo- 
cal experimental techniques we extend the algebraic re- 
sults for the single extended impurity to a finite im- 
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purity concentration m within the self-consistent T- 
matrix (SCT) approximation. In the SCT approximation 
we obtain the impurity position averaged self-energy as 
t(k,u) = n i f kM (Lj). With G(k,k,w) = [(Go(k^))- 1 - 
E(k, uj))' 1 replacing Go in the above exact formulas for 
the single impurity T-matrix we solve the resulting self- 
consistency equations for G by iteration without the 
usual assumption of particle-hole symmetry. Note that 
£ = X^=oi3^^' i- e - t nc T\ component of the self- 
energy is finite for a > 0. 

Fig. 4 shows the expected result for the SCT DOS for 
a low impurity concentration — 1%. At first sight the 
zero frequency DOS increase with turning on a might 
seem weak. However, within the chosen parametrization 
a = 0.5 translates into Vi/Vo — 1/16; i.e. if the nn 
impurity potential strength is only 6% of the impurity 
center value and accompanied by a suppression of the 
hopping amplitude to the impurity center, N(u> = 0) 
increases by ~ 15% which is indeed significant. 

In summary, we have solved a physically motivated 
extended impurity potential model in a d x 2_ y i supercon- 
ductor. The finite range of the potential is argued to 
be dynamically generated around a local impurity from 
strong electronic correlations in the host system and thus 
relevant for high-T c superconducting materials. While 
the model itself describes a static potential with a se- 
lected parameter relation it has the advantage to allow 
for an algebraic solution by decomposition into s— , p-, 
and d-wave scattering channels. The significant enhance- 
ments of the impurity induced resonances by the finite 
potential range clearly improve the proposed possibility 
for probing the highly anisotropic local DOS structure 
near the impurity by STM experiments and may thus 
serve as an alternative fingerprint of the gap symmetry 
in cupratc and other unconventional superconductors. 
Furthermore, due to the physical relevance for the de- 
scription of impurity effects in high-T c superconductors 
some of the previous model studies with s-wave scat- 
tering from local impurity potentials may be reconsid- 
ered for the effects of an extended potential range for 
which the presented model can serve as an algebraically 
tractable tool. 
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FIG. 4. Impurity position averaged DOS from the self- 
consistent T-matrix approximation for an impurity concen- 
tration of Hi = 1%. For the meaning of the different lines see 
Fig.2. 
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